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The main characteristics of the linear Darrieus-Landau instability in the laser ablation flow are investigated. 
The dispersion relation of the instability is found numerically as a solution to an eigenvalue stability problem, 
taking into account the continuous structure of the flow. The results are compared to the classical Darrieus- 
Landau instability of a usual slow flame. The difference between the two cases is due to the specific features of 
laser ablation: high plasma compression and strong temperature dependence of electron thermal conduction. It 
is demonstrated that the Darrieus-Landau instability in laser ablation is much stronger than in the classical case. 
In particular, the maximum growth rate in the case of laser ablation is about three times larger than that for slow 
flames. The characteristic length scale of the Darrieus-Landau instability in the ablation flow is comparable 
to the total distance from the ablation zone to the critical zone of laser light absorption. The possibility of 
experimental observations of the Darrieus-Landau instability in laser ablation is discussed. 



I. INTRODUCTION 



Inertial confinement fusion (ICF) is believed to be one of the promising energy sources in the 21 st century. The aim of ICF 
is to compress a plasma target to densities and temperatures high enough to trigger a thermonuclear reaction. Despite a great 
technical development and progress in power supplies during the last decades fljj] , hydrodynamic instabilities remain the limiting 
factor in fusion performance and efficiency. In this respectjthe most difficult obstacle in achieving ICF is the Rayleigh-Taylor 
(RT) instability, which arises because of target acceleration SB SHE 8, 9]. Still, the RT instability is not the only instability 
of importance in ICF; for example, the laser generated plasma is also subject to the so called Darrieus-Landau (DL) instability 
lIlJSMM, El El El El El- Traditionally, the DL instability is known as the hydrodynamic instability of a slow flame 
118, 19, 20l l2lll . which makes the flame front corrugated and increases the burning rate. Flame is the most typical example of a 
deflagration wave, which is a front propagating due to energy release and thermal conduction. The ablation flow generated by 
laser radiation on a plasma target is also a deflagration wave H [H |23| . Other interesting examples of deflagration come from 
astrophysical applications like the big bang model and type I supernovae ll2^l2^l2^l27il28il29tl . On the basis of similar physical 
properties of deflagrations, one should expect the DL instability to develo p in laser ablation. There has been much interest in the 
DL instability in ablation flows, see e.g. Refs. ll5l ITol ITTl IT2I Il3l IT41 1 1 5L 1 1*61 IT7h . However, the most important features of the 
instability has remained unclear up until now. The purpose of the present theoretical work is to answer some of the key questions 
concerning the DL instability in laser ablation and to indicate conditions for experimental observation of the instability. 

The classical theory of the DL instability considers an incompressible flow generated by a slow flame lfl8l El Hi Eh. 
However, the parameters of laser ablation flows in ICF are markedly different from those of slow combustion. One of the most 
distinctive features of deflagration in ICF is that the plasma velocity reaches the isothermal sound speed at the critical surface 
of laser light absorption. For this reason, laser ablation corresponds to the so-called Chapman- Jouguet (CJ) deflagration, which 
is the fastest propagation regime possible for a deflagration front. Thus, in order to describe the DL instability in laser ablation 
one has to take into account strong compression of the plasma flow with relatively high local values of the Mach number. 
Another important property inherent to the plasma flow is the strong temperature dependence of electron thermal conduction. 
Electron thermal conduction determines the total thickness and the internal structure of laser deflagration. It is expected that 
the strong plasma compression and the electron thermal conduction properties will significantly influence the properties of the 
DL instability in a laser generated plasma, making it markedly different from the case of slow flames. A number of papers was 
devoted to the linear and nonlinear stages of the DL instability in ICF d [U [H EH EH EH EHl- Still, these papers did not 
answer the most important questions concerning the DL instability in laser ablation. In particular, we would here like to outline 
the following three questions within the problem: 



1. How strong is the DL instability in laser ablation in comparison to the classical case, i.e., is it stronger or weaker than the 
DL instability developing at a slow flame front? 

2. What is the characteristic length scale of the instability development? This question is especially important from the 
experimental point of view, since the answer determines the target size, for which the DL instability may be observed. 

3. What is the outcome of the DL instability at the nonlinear stage? 

There have been numerous attempts to answer the first question within the model of a discontinuous deflagration front in a 
compressible gas/plasma flow ITol [Til [l6l l30l [3lll . Unfortunately, the discontinuous model encounters the deficit of matching 
conditions at the deflagration front; the number of unknown variables exceeds the number of conservation laws at the front by 
one. The problem was encountered first within the studies of the RT instability in laser ablation, see Refs. 13t|5i|9J]; more detailed 
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discussion on this subject may be found in Ref. Il6ll . In the classical incompressible limit of the DL instability at a slow flame 
front, the extra condition is given by the so-called DL condition of a constant deflagration speed with respect to the cold gas. This 
condition may be proven rigorously, see Ref. lEol l2ltl . A counterpart of the DL condition for a compressible flow is not obvious 
and the solution to the problem turned out to be sensitive to this assumed extra condition, so that the different analytical theories 
proposed led to qualitatively different results. For this reason, it was unclear if strong plasma compression in laser ablation made 
the DL instability stronger, or weaker, or only produces minor changes in the instability strength. Within this context, Refs. 
I32l,l33ll deserve special attention, since these papers considered influence of gas compression on the DL instability taking into 
account a continuous structure of the deflagration front. Though Refs. 11321 13311 were devoted to normal combustion, and the 
respective results cannot be extrapolated directly to the DL instability in laser ablation, the method used in these papers may 
still be used for the ablation studies. This method eliminates the deficit of boundary conditions and allows for investigating the 
properties of the DL instability in laser ablation. Here we employ the methods of Refs. [32, 33] to study the linear stage of the 
DL instability in a laser plasma, thus answering questions 1 and 2 of those outlined above. 

In the present paper we investigate main characteristics of the linear DL instability in the laser ablation flow. We find the 
dispersion relation of the instability numerically as a solution to an eigenvalue stability problem taking into account the contin- 
uous structure of the flow. We compare the results to the classical DL instability of a usual slow flame. We show that difference 
between the two cases is due to two specific features of laser ablation: a high plasma compression and a strong temperature 
dependence of the electron thermal conduction. We further demonstrate that the DL instability in laser ablation is much stronger 
than in the classical case. In particular, the maximal growth rate of perturbations in laser ablation is about three times larger than 
for slow flames. The characteristic length scale of the DL instability in the ablation flow is comparable to the total distance from 
the ablation zone to the critical zone of laser light absorption. We discuss the possibility of experimental observations of the DL 
instability in laser ablation. 

II. BASIC EQUATIONS AND THE STATIONARY SOLUTION 

We describe the ablation plasma flow using hydrodynamic equations of mass, momentum and energy conservation 

^+V.(pu)=0, (1) 

p^ + (pu-V)u + VP = 0, (2) 
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and the equation of state of an ideal gas 
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P= L C P pT, (4) 
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where 7 = 5 /3 is the adiabatic exponent, Cp and Cy are the heat capacities at constant pressure and volume respectively. Electron 
thermal conduction K depends on temperature as K = K C (T /T c ) 5 / 2 , where label "c" refers to the critical surface of laser light 
absorption. Laser light absorption brings energy into the plasma and, together with thermal conduction, it drives the flow. 
Absorption takes place when plasma frequency is equal to the laser frequency 
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which determines plasma density p c at the critical surface (here M is plasma mass per one electron). Decrease of the laser light 
intensity due to absorption may be described as |Q] 

f=Kl, (6) 
dz 

with the absorption coefficient 
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The absorption coefficient diverges at the critical surface, and, therefore, the process of energy absorption is strongly localized at 
the surface. In the studies of ablation flow and the RT and DL instabilities in the flow, the energy release is typically presented by 
5-function JH HI [HI [HI HH . Such replacement is possible since the instabilities develop on the length scales much larger than 
the region of energy release and involve bending of this region as a whole without changing its internal structure. In the present 
paper we solve the problem of the DL instability numerically. In the numerical solution, it is more convenient to imitate the 
5-function by a transitional zone of finite width determined by some continuous function Q.r of energy gain in the flow included 
into Eq. (0. Here, we chose the function in the form suggested in Ref. II 1611 



The function Q.r was constructed taken into account similarity with the Arrhenius law in combustion, where j3 plays the role 
of the scaled activation energy and n is similar to the reaction order. It is well-known that the Arrhenius law provides strong 
localization of energy release for any reasonable n in the case of sufficiently large J3 >> 1, see lfl9ll20h . In combustion science 
the Arrhenius reaction is sensitive to temperature changes [ 19, 20]. Here we construct the function of energy gain Q.r sensitive to 
density variations, as it takes place in laser ablation. Choosing large parameter j3 — > °° we obtain energy gain strongly localized 
at the critical surface with p — * p c . On the other hand, the numerical solution demands a finite width of the zone of energy 
release and the finite value of the parameter p\ In most of our calculations we use j3 = 90, n = 2, We also investigate sensitivity 
of the physical results to the choice of these parameters. We demonstrate that these parameters have minor influence upon the 
properties of the DL instability and at high values of j3 this influence vanishes. The function Q,r given by Eq. (O allows a planar 
stationary solution consisting of two uniform flows of cold heavy plasma (label "a") and hot light plasma (label "c") separated 
by a transitional region, which is the deflagration front. The labels "a" and "c" originate from the ablation and critical surfaces 
in the laser deflagration. Typical internal structure of the deflagration front is illustrated in Fig. 1 . The described geometry is 
common in the theoretical studies of the RT and DL instabilities of the ablation flow , though it does 

not take into account the rarefaction wave in the hot light plasma beyond the critical surface. As the main advantage of such a 
choice, we may consider different values of the Mach number in the light plasma, which would be impossible with a rarefaction 
wave. Changing the Mach number, we can go over continuously from the case of classical incompressible DL instability of a 
usual flame front to the case of laser ablation with strong influence of plasma compression. We stress that the main purpose of 
the present paper is to compare the DL instability in laser ablation to the classical case. Therefore, in general, we will discuss 
stability of a deflagration front, which covers the cases of usual flames (slow combustion) and laser ablation as two asymptotic 
limits of negligible compression effects and ultimately strong plasma compression. 

We start our analysis with describing internal structure of the stationary planar deflagration (laser ablation) flow; this is the first 
step in the stability analysis. Figure 1 illustrates plasma density, temperature and energy release in the flow obtained numerically 
as described below. The deflagration front in Fig. 1 propagates to the left with constant velocity U a (the ablation velocity) in 
the negative direction of z-axis. Velocity of a usual flame is determined by the rate of energy release and thermal conduction. 
Ablation velocity is determined by the critical density and the laser light intensity 02211 . We adopt the reference frame of the 
deflagration front. In that case the front is at rest, but the cold heavy plasma flows to the right with uniform velocity u z = U a , 
undergoes transition in density and temperature in the deflagration wave, and, finally, the hot light plasma gets drifted away with 
uniform velocity u z = U c . 

Equations of mass and momentum transfer ([D, (f2]l may be integrated for the planar stationary deflagration as 



Q. R = Q.(p-p c ) n exp( 



(8) 



pu z = p a U a = p c - 
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P + puj = P a + PaPl = Pc + pcUc 2 . 



(10) 



One of the main dimensionless parameters in the problem is the expansion factor 







Pa 
Pc 



Hi 



(ID 



which shows density drop from the original cold plasma to the critical surface. Both in flames and laser ablation, the expansion 
factor is rather large, © = 5 — 10, see J3H3l. In the case of ablation flow, the laser light frequency determines the critical density 
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and the expansion factor. The other important parameter is the Mach number in the light plasma (gas) corresponding to the 
adiabatic sound 




The Mach number is negligible in the classical case of usual flames, and it may be taken zero with a very good accuracy. On the 
contrary, laser ablation provides an ultimately large value of the Mach number possible for a deflagration flow. In laser ablation, 
the isothermal Mach number is equal unity p c U 2 /P c = 1 in the light plasma, and we have Ma 2 = 1/7 for the adiabatic Mach 
number. In the present work we consider a general case of a deflagration front with an arbitrary Mach number changing within 
the limits < Ma 2 < 1/y. The internal structure of the deflagration front follows from the stationary equation of energy transfer 
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Characteristic width of the front is determined by thermal conduction in the hot region 



L c = --£-. (14) 
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The problem involves one more parameter of length dimension 
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related to thermal conduction in the cold flow. Because of the strong temperature dependence of electron thermal conduction, 
these two length scales are quite different L c /L a — (T c /T a ) 5 / 2 . For example, for the temperatures ratio T c /T a = 6, the length 
scales differ by two orders of magnitude L c /L a ss 88. The strong difference in these length scales is one of the specific features 
of laser ablation in comparison with usual flames. 

We introduce dimensionless variables for plasma density, temperature, velocity and coordinate 
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Then we can rewrite Eq. ( fT3T > as 
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where A = L 2 Q.p? {KcT c )~ l is an eigenvalue of the stationary problem. The relation between temperature and density in a 
deflagration flow follows from Eqs. (0, ( fTOb 



In the incompressible limit of usual flames, Ma 2 << 1, this relation is reduced simply to (p9 = 1, so that temperature ratio T c /T a 
is determined by the expansion ratio, T c /T a = 0. In the case of strong gas compression, these two values differ considerably. 
For example, in the case of laser ablation with the critical Mach number Ma 2 — 1 /y, we find from Eq. (fT8l that 



T c ® 2 

— = • (19) 

T a 20-1 v ; 

For high values of density drop, 2© >> 1, temperature ratio is about twice smaller than the expansion factor, T c /T a w 0/2. 
Because of the reduced temperature ratio, the effect of two different length scales Eqs. ( TBI . ( fT5l ) in the ablation flow is expected 
to become weaker than in a similar incompressible flow. With temperature ratio changing from T c /T a = in the incompressible 
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case to T c /T a w 0/2 in the ablation flow we find the ratio of length scales L c /L a — (T c /T a ) 5 / 2 decreasing by the factor of 
©V 2 w 5.7. For this reason, the profiles of density and temperature are expected to be much smoother in the ablation flow in 
comparison with an incompressible counterpart. 

We solve Eq.dTTb together with Eq. ( fT8l l numerically for different values of the Mach number. Typical solution to the problem 
for incompressible case is shown in Fig. 1. First of all, we have to make sure, that the ablation front structure is not sensitive to 
our choice of the energy gain function Q.r . We investigate dependence of the density and temperature profiles on the parameter 
/3. Figure 2 shows density and energy release for = 6, n = 2, Ma c = 0.5, j3 = 20; 60; 140. Density profiles for these three 
values of /3 coincide almost everywhere except for the zone of energy gain. Still, even inside this zone, the difference between 
the density profiles is minimal, and cannot be recognized on Fig. 2. Similar result holds for the temperature profiles. Thus, 
we can conclude that parameter j3 does not affect the density and temperature profiles in the ablation front. At the same time, 
parameter j3 influences strongly the profile of energy release itself: the respective plots become much more localized with /3 
increasing. Thickness of the critical surface should be zero in the hydrodynamic description of the ablation flow. Therefore, 
finite width of the zone of energy gain implies certain inaccuracy of the numerical solution. The level of inaccuracy of the model 
may be evaluated as the characteristic width of the energy gain zone in the dimensionless variables (scaled by the total thickness 
of the deflagration front L c ). As we can see in Fig. 2, half width of the energy peak is rather wide for j3 = 20; it takes more than 
20% of the whole deflagration front. As we increase j3, this width becomes smaller and for j3 > 90 the inaccuracy is less than 
5% . The other parameter of the energy gain function in Eq. (JHJ, «, affects the shape of the energy release peak only slightly 
even for moderate values of j3. For j3 = 90 the parameter n has a negligible effect on the physical results. Thus, Eq. ([8]l may 
imitate the energy gain in the ablation flow quite well; in all following numerical solutions and figures j3 = 90 and n = 2. 

As the next step, we solve Eq. ( fTTb for different values of the Mach number. Figure 3 shows profiles of density and energy 
release, for Ma c = 0; 0.5; 0.75, other parameters being fixed as © = 6, j3 = 90, n = 2; Figure 4 presents the respective temperature 
profiles. In the case of incompressible flow, Ma c — 0, the density profile demonstrates clearly the effect of two length scales, Eqs. 
(O, (fTBT l, produced by the temperature-dependent thermal conduction. The profile is rather smooth in the hot region close to the 
critical surface, and it becomes sharp in the cold plasma close to the ablation surface with large density gradient. In fact, it is this 
property, which allows distinguishing two effective surfaces in the flow: the ablation and critical surfaces JHfSHl- In the case 
of strong compression with Ma c = 0.75, the density profile becomes much smoother. Smoothing of the density profile concerns 
the region of cold plasma mainly. Another specific feature of the density profile at high values of the Mach number is shown at 
the insert of Fig. 3. With the Mach number approaching the maximal possible value Ma 2 . = 1/y, we observe development of 
another mini-region of relatively high density gradient close to the critical surface. This effect may be also obtained analytically 
from Eq. ( TT8l >. Close to the critical surface, expanding density and temperature in power series with respect to <p — 1 << 1, 
1 — << 1, we obtain the relation ip-ls; ^1 — 6, or p/p c ■— 1 ~ J 1 —T/T c in the dimensional values. Taking energy gain 
in the form of 5(z)-function, we have temperature achieving the final value T c at finite point z = smoothly. This leads to the 
square-root singularity in the density gradient. We can observe the trace of such a singularity in the insert of Fig. 3, though 
smoothed because of the finite width of the energy gain zone. We also observe changing shape of the energy gain with Mach 
number. The energy gain zone becomes much thinner at high values of the Mach number. This happens because the energy 
release Eq. ([8]) is sensitive to the density profiles. Sharp gradients of density at Ma 2 = 1/y make the zone of energy gain sharper 
as well. As a result, the model Eq. dS) works much better at high values of the Mach number corresponding to the laser ablation 
flow. Figure 4 demonstrates that temperature profiles become also smoother with increasing Mach number, which is similar to 
density profiles. Besides, temperature ratio at the ablation and critical surfaces decreases with increasing the Mach number, as 
we demonstrated in Eqs. ( fT8l l, ( fT9l ). The numerical solution for the planar stationary flow illustrated in Figs. 1-4 provides the 
basis for the stability analysis performed in the next section. 

III. LINEARIZED EQUATIONS 

We solve the stability problem for small perturbations of any value (j) in the form: 

<j){x,z,t) = <j>(z) + ij>(z)exp(ot + ikx), (20) 

where the first term in the right-hand side of < f20b stands for the stationary flow, the second term describes linear perturbations, 
a is the instability growth rate and k — 2n/X is the perturbation wave number. In general, a may have both a real part (growth 
rate) and an imaginary part (frequency). However, in the case of the DL instability a is only real; the instability develops when 
(7 is positive. 

The linearized system (1-3) takes the form 
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Similar to 15[], we introduce the dimensionless perturbations of mass flow j, transverse velocity v, temperature 9 and dynamic 
pressure II as 
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with the scaled wave number and the perturbation growth rate 
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Then the linearized system (f2TT>-(f24b is 
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where w = 9 — yMa 2 u and the coefficients in Eq. d30b are 
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As the boundary conditions to the system d27li-(l30il. we demand that perturbations vanish at infinity in the uniform flows of cold 
plasma ahead of the ablation zone and the hot plasma behind the critical zone of energy gain. The coefficients in Eqs. (l27Ti-(l30t 
are constant in the uniform flows. This allows us writing down the perturbations in an exponential form = 0exp(jU<i;), 
where jj. > in the heavy plasma (£, — > — °o) and jj, < in the light plasma (£, — > °o). In some particular simplified limits the 
structure of the perturbation modes in the uniform flows may be written analytically, e.g. see fl6|]. However, in the present case 
we can do it only numerically, which becomes another step in the numerical solution to the problem. 



IV. THE ALGORITHM FOR THE NUMERICAL SOLUTION 



The derived system (I27li-(l30l> is rather complicated and it may be solved only numerically. We introduce an auxiliary variable 
\f/ = 5 / 2 d9/ d% to obtain two differential equations of the first order instead of Eq. d30b as 

d9 5 1 d0 ~ , /9 ~ 

d-r~2-e^ e+0 V ' (36) 
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Thus we have a system of five differential equations of the first order with the scaled growth rate S as an eigenvalue. The purpose 
of the solution is to find the dispersion relation S = S (K). The system may be written in a matrix form 
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F = 



w 


-K<p 


w 


W 







W 


-S(p 


W 


W 







-s 


-K 











(39) 











5 1 de 

2 c/§ 


0-V2 




A) 


A v 


An 


A e 


Ay/- 





The numerical solution consists of the following steps. First we look for the density profile in the stationary solution to Eq. 
Equation (T% is integrated numerically from the uniform flows of cold and hot plasma to the central parts of the deflagration 
transitional region. Both solutions are matched at a certain density value; we checked that the physical results do not depend on 
the choice of the matching point. When density distribution is known, temperature and velocity profiles are determined using 
Eqs. © and (1181 . The respective numerical problem involves several length scales, which creates an additional difficulty in 
the solution. In the case of highly compressible flow the peak of energy release is extremely sharp; but we have to resolve it 
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properly, since it is used afterwards in the solution to the stability problem. As the next step, we find modes 0(4) = exp (ju4) 
in the uniform flows, which determine the boundary conditions for Eqs. ( f27l > - ( f29b . In the uniform flows, Eq. ( l38l is a system 
of linear ordinary equations with constant coefficients. In order to find the factors /I, we solve 

|F-Eji|=0, (40) 

where E is the unit matrix. Thus we obtain an equation for ji in the form of a polynomial of the fifth order. In the incompressible 
case this equation may be solved analytically. Five different roots correspond to the vorticity mode, two sound modes and two 
modes of thermal conduction and/or energy gain 15lll6|,|34]]. Solving Eq. (l40l > numerically we find five modes with two positive 
and three negative values taking ji > for 4 — * —°° and < for 4 — > °°- Finally, we integrate the system d3~8l numerically 
in the transitional region of the deflagration flow. We perform the numerical integration two times from the right-hind side and 
three times from the left-hand side with boundary conditions determined by different modes. We match these five solutions at a 
certain point between the maximum of the energy release and the ablation zone of sharp density gradients. Again, the physical 
results do not depend on the choice of the matching point. Then we obtain a matrix consisting of twenty five values describing 
flow perturbations for five modes. Taking the determinant of this matrix equal zero, we find the dispersion relations' = S (K). 

V. RESULTS AND DISCUSSION 

Before presenting the numerical solution to the stability problem, we explain the physical results we are looking for. In the 
classical case of an infinitely thin flame front propagating in an incompressible flow, the instability growth rate was obtained by 
Darrieus and Landau lfl8l[T9ll as 

a = rU a k, (41) 
where the coefficient T depends on the expansion factor © only 

® V®+1- ■ ( 42 ) 
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Traditionally the DL instability growth rate is scaled by the front velocity U a instead of U c — &U a used as the velocity unit in 
Sees. II, III. Here we follow the tradition and use U a for scaling when presenting the results. The approach of an infinitely 
thin front of the classical solution holds for long wavelength perturbations kL c < < 1 . Taking into account gas compression, 
we should expect the dispersion relation for an infinitely thin deflagration front in the same form as Eq. (l4TT i. but with the 
coefficient T depending on the Mach number, T = T(&,Ma c ). More general solution to the stability problem takes into account 
finite thickness of the deflagration front. Finite thickness of the deflagration front leads to stabilization of the DL instability at 
sufficiently short wavelengths. In the case of usual slow flames of finite thickness, the analytical solution to the problem may be 
found e.g. in l2lll . Written in the form of Taylor expansion in relatively small perturbation wave number, kL c << 1, the solution 
may be presented as 

a = TU a k(l-k/k cut ), (43) 

where k cut is the cut-off wave number, X cm = 2n/k cut is the cut-off wavelength. The approximation of the small cut-off wave 
number k cut L c < < 1 holds with reasonable accuracy for usual flames, e.g. see the review [20] . The cut-off wavelength is propor- 
tional to the thickness of the deflagration front, X a , t « L c . In the incompressible flow, the coefficient of proportionality depends 
on the expansion factor © and on the type of thermal conduction. Particularly, in the case of K °= T 5 / 2 and an incompressible 
flow, the theory [21] predicts 



Km 47T0 



Lc ©-1 



0+1 



(44) 



For typical expansion factors = 6 — 8, Eq. d44b predicts the cut-off wavelength X cut « 6L C . For comparison, the DL cut- 
off for usual flames is considerably larger, being about k cut ~ 20L C , see [20]. According to the dispersion relation ( |43l , there 
is a maximum of the instability growth rate (J max achieved at a certain finite perturbation wavelength A max . Equation d43l 
predicts the wavelength of the maximum to be twice larger than the cut-off wavelength, A max = 2X cut . Taking into account gas 
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compression, one should expect that all these values depend on the Mach number, Ma c . Therefore, the purpose of the present 
work is to investigate dependence of the parameters Y, X a ,t, A max , cr max on the Mach number Ma c changing within the interval 
< Ma c < 1 / y 7 /. We demonstrate below that this dependence is quite strong. 

Figure 5 presents our numerical solution to Eqs. d27ll-(l30ll. that is the scaled instability growth rate versus the perturbation 
wave number. The dispersion relation is shown for different values of the Mach number, Ma c = 0; 0.5; 0.65; 0.7; 0.73, with the 
expansion factor = 6 and a strongly localized energy release for /3 = 90, n = 2, see Sec. II. As we can see, the DL instability 
becomes much stronger with increasing the Mach number; plasma (gas) compression provides a strong destabilizing effect. This 
result looks very similar to one obtained for a flame in a compressible flame in Refs. IB21I33I1 . The destabilization concerns all 
parameters of the instability: the factor F in Eq.fiTI). the maximal instability growth rate, <7 ma x, and the cut-off wave number, 
X cut . These results are presented in Figs. 6-8, respectively. Figure 6 shows the factor T, which determines strength of the DL 
instability in the case of an infinitely thin front. The solution is obtained for the expansion factors © = 6, 10. According to Fig. 
6, the DL instability for an infinitely thin ablation front (with ultimately strong plasma compression) is almost twice stronger 
than in the incompressible flow. The ratio of T-factors obtained for Ma c = 0.73 and Ma c = is about 1.8 for © = 6, 10. Some 
analytical theories for the DL instability in compressible flows predicted also increase of the growth rate with the Mach number 
fl3C ,13111 . Still, the numerical results show a noticeably stronger increase than any analytical theory proposed so far. As explained 
in itlol . the stability problem of a discontinuous deflagration front encounters a deficit of matching conditions at the front: the 
number of unknown values exceeds the number of matching equations by one. In order to overcome the obstacle, different 
additional matching conditions were suggested in different papers ITol [Til [l6l 1351 13111 , Solution to the problem turned to be 
quite sensitive to the choice of the extra condition; but the extra conditions suggested were merely assumptions. This trouble 
does not happen in the numerical solution, since the numerical solution considers continuous transition from heavy cold plasma 
to light hot one. For this reason, the numerical solution provides also a test for the suggested analytical solutions. In this work 
we are not going to criticize the previous analytical theories. Instead, we will try to extract the best ideas suggested so far in the 
theoretical papers ITcil [Til [T6l l30l [3lll to obtain the analytical solution reasonably close to the numerical one. For that purpose 
we take the basic elements in the analysis tig] and complement them by an extra matching condition identical to that of the DL 
theory of the incompressible flow |[T8ll . In the dimensional variables the DL matching condition is 
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where / is perturbation of the discontinuous front position. The same type of matching condition was assumed for the DL 
instability in a compressible deflagration/ablation flow in fTol [TTl l30ll . Reproducing calculations of lfT6Tl with the extra condition 
d45l l we obtain the following equation for F 

1 - pe-i ml r„, -1 a + effi* _ „ 

1 - ©Ma 2 , F- TJ e ©r 1 — ©Afa 2 

where 



77„ = ^1+Mfl2(@2r2-1), Tj c = ^l+M fl 2(r 2 -1) (47) 
and the Mach number in the cold plasma 

9 Ma}. 

Maf, = ; — —s t ■ (48) 

" © + /(©- l)Ma 2 V ' 

Influence of the Mach number may be illustrated in the limit of small plasma compression, Ma << 1, using Taylor expansion of 
Eq. d46i >. In that case 

r = r (i + /3Mfl c 2 ), (49) 

where Fq corresponds to DL solution for incompressible case and j3 is determined as 

2©(r (© + 2) + i) >Q 
(© + i) 2 (r (© + 1) + ©) 
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Positive factor j3 indicates increase of the instability growth rate with plasma compression. The numerical solution to Eq. (146b 
together with Eq. d49l are presented in Fig. 6 by the solid and dashed lines respectively. These lines do not provide a perfect 
agreement with the numerical results; still the difference between the theory and the numerical solution is acceptable, about 15% 
. For comparison, Fig. 6 shows also the instability growth rate predicted in Q30I1 . by dash-dotted lines. The theory in [30] differs 
much stronger from the numerical solution, approximately by 35% . 

Figure 7 shows the maximal instability growth rate versus the Mach number for the expansion factors = 6, 10. As we can 
see, in the case of laser ablation with Ma c = 1 / ^fy the maximal instability growth rate is about three times larger than in the 
incompressible case. These results agree well with the previous numerical calculations of Ref. j32tl for flames in a compressible 
flow. It is interesting that the maximal growth rate shows only minor dependence on the Mach number for a rather wide range of 
this parameter, Ma c < 0.5. The strong dependence of (7 max on Ma c takes place only when the Mach number starts approaching the 
limiting value Ma c = 1 /y 7 / inherent to laser ablation. We observe a similar tendency in Fig. 8, which presents the cut-off wave 
number versus the Mach number. Again, the cut-off wave number is about twice larger for laser ablation in comparison with the 
incompressible case of zero Mach number. Figure 9 compares the cut-off wavelength found numerically for different values of 
the Mach number to the theoretical prediction Eq. (144-b . We remind that Eq. (l44l) follows from the theory [21j in the limit of 
an incompressible flow and thermal conduction depending on temperature as K T 5 I 2 . As we can see, the theory d44l > provides 
a reasonable prediction for the cut-off wavelength in the case of zero Mach number; the difference between the theory and the 
numerical solution is about (15-25)% . As we increase the Mach number, the DL instability becomes stronger and the cut-off 
wavelength decreases considerably. For example, taking the expansion ratio © = 7 we find the cut-off wavelength X cut sa 2AL C 
for laser ablation, X cut w A.%L C for the incompressible case of Ma c — and X cut 6L C predicted by the analytical formula d44b . 
Thus, the numerical solution predicts the DL instability in laser ablation on the length scales larger by the factor of 2.4 than the 
distance from the ablation zone to the critical zone. These length scales are very small when compared to the respective ratio 
k C ut/L c ~ 20 for usual flames. However, these length scales are extremely large in comparison with the length scales typical 
for the RT instability in inertial confined fusion fH |(| HI] • For this reason, in order to observe the DL instability experimentally 
one has to take special precautions eliminating the RT instability in the flow. One of the possible options is to use a sufficiently 
large target of the thickness exceeding the distance L c at least by an order of magnitude. In that case target acceleration is minor, 
which leads to the negligible RT instability, while the DL instability has sufficient space to develop. In addition, one has to 
remember that experimental observations typically concern the fastest growing perturbations of the wavelength A max and larger, 
not the cut-off wavelength X cut . In the case of usual flames these two length scales are related as Amax ~ 2Ac«/ ) since the whole 
dispersion relation may be described with a good accuracy by two first terms in Taylor expansion in kL c << 1, see Eq. d43l . 
The ratio A, nax /A CU f becomes somewhat different for the DL instability in laser ablation. Figure 10 shows the ratio Amax/A c „; 
obtained numerically for different values of the Mach number, and compares it to the classical case of Amax/A c „ r = 2. We can see 
that the wavelength corresponding to the maximal instability growth rate is about Am ax ~ 1 -8A C „ ; in laser ablation. The deviation 
indicates that next order terms become important in the Taylor expansion of the instability growth rate in kL c << 1. Still, the 
deviation is not too large. 

Finally, we have to check that our model for the energy gain in the deflagration/laser ablation flow does not influence the 
physical results obtained. In order to validate the model we investigated influence of the parameters j3 and n of the energy 
gain on the DL dispersion relation. Numerical calculations for j3 = 90 and different values of n show negligible variations of 
all parameters of the instability: the T-factor, the maximal instability growth rate and the cut-off wavelength. We also took 
n = 2 and varied /3 within the limits between 20 and 140. For example, taking Ma c — 0.65 we find the maximal growth rate 
Cmax^c/t/u = 0.74 for /3 = 20, OmiaLc/Ua = 0.84 for J3 = 90 and a ma _ x L c /U a = 0.86 for j3 = 140. These calculations indicate 
that the continuous numerical model for the energy gain function £2r in Eq. <[8J brings inaccuracy of only few per cent, about 
3% , into the numerical solution for j3 = 90 used in the present paper. Investigation of the cut-off wavelength for different j3 and 
n leads to similar conclusions. 



VI. CONCLUSIONS 

In the present paper we investigated the DL instability in an ablation flow and compared the results to the classical case 
of a slow flame. Unlike the normal flame, laser ablation is characterized by the strongest plasma compression possible for 
a deflagration wave. Another specific feature of laser ablation is the strong dependence of electron thermal conduction on 
temperature. We demonstrate that the DL instability in laser ablation is much stronger than in the classical case. In particular, 
the maximal growth rate in the ablation flow is about three times larger than in the incompressible case. Moreover, the cut-off 
wavelength changes drastically as we go from the classical case of an incompressible flow to the ablation flow. The cut-off 
wavelength is also strongly influenced by the temperature dependence of thermal conduction. It is known that the DL instability 
for usual flames develops on quite large length scales exceeding the flame thickness by almost two orders of magnitude [20]. 
In contrast to this, the characteristic length scale of the DL instability in the ablation flow (e.g., the cut-off wavelength) is 
comparable to the total distance from the ablation zone to the critical zone of laser light absorption, L c . Still, even these values 
are large from the point of view of possible experimental observations of the DL instability in laser ablation. We note that the RT 
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instability in inertial confined fusion develops on length scales much smaller than L c . For this reason, the DL instability may be 
observed only if the accompanied RT instability is suppressed. This may be achieved, for example, for sufficiently large targets 
of thickness much larger than the distance L c from the critical to the ablation zone. 
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FIG. 6: The T factor versus the Mach number squared for © = 6 (circles) 0=10 (squares). The solid lines depict Eq. i46l . the dashed lines 
corresponds to Eq. i49\ . the dash-dotted lines present the result of Ref. 13011 . 




FIG. 8: The cut off wave number versus the Mach number, © = 6; 10. 
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FIG. 9: The cut off wavelength versus the expansion factor forMa e = 0; 0.65; 0.73. The dashed line shows the analytical formula Eq. l !44t 
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